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Abstract 

Recently Sun and Gong proposed a new constant-pressure molecular dynamics method for finite 
systems. In this paper, we discuss the current understanding of this method and its technique 
details. We also review the recent theoretical advances of nano-system under pressure by using 
this method. 



/. Introduction 



Nowadays, the molecular-dynamics (MD) simulation, widely used in chemistry, physics, 
and materials sciences, is considered as a standard and powerful tool for investigating the 
structures and properties of matters in atomic scale.- As an important improvement, the 
constant-pressure MD(CPMD) proposed by Andersen,- and subsequently extended by Par- 
rinello and Rahman,- has opened a crucial window to explore systems under the pressure 
and tensions. Over the past decades, CPMD plays a key role for our understanding of many 
phenomena relevant to high pressure experiments in atomic scales. 

Although the traditional CPMD has archived great success for bulk matters at high- 
pressure conditions, it fails to be directly used for finite size systems and for systems with- 
out regular shapes, such as nanocrystals, where the boundary is hard to describe. Mo- 
tivated by the experimental work for the molecular, low- dimensional, biological system^ 
and nanocrsystals under pressure,- three theoretical groups proposed the CPMD for fi- 
nite systems. To achieve the goal, Martonak, Molteni and Parrinello have made the first 
step by directly extending the traditional CPMD.- Hereafter we call this method as the 
directed method. This method directly mimics a real high-pressure experiment, i.e. the 
system keeping at constant pressure through exchange of linear momentum with environ- 
ments (in their paper the environment is the pressure transmitting liquid). According to 
Martonak, Molteni and Parrinello,- a target cluster is immersed into a well-chosen pressure- 
transmitting liquid, the whole system (liquid+cluster) is simulated using Parrinello- Rahman 
CPMD. The Parrinello-Rahman CPMD is a well developed technique, so the key point of 
using this method is the choice of pressure-transmitting liquid, specifically the interaction 
between liquid-liquid and liquid-cluster. In the real application, these interactions should 
be set to prevent the liquid from being inside the cluster, and from phase transition hap- 
pening during the simulation. Furthermore one should pay more attentions to the choice of 
pressure-transmitting liquid when the target cluster being in liquid state due to the diffusiv- 
ity In their paper, Martonak, Molteni and Parrinello have used a classical repulsive liquid, 
and with classical interactions between the cluster and the liquid.- The main drawback of 
this method is the number of the pressure-transmitting liquid atoms should much larger 
than that of the cluster. Additionally, the direct method is suffered the same problem as 
the original CPMD by the artificial mass associated with the piston. 
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Kohanoff, Caro and Finnis presented another method by introducing the stochastic Brow- 
nian forces to each surface atoms (hereafter stochastic method).- This method can be con- 
sidered as simplified version of the direct method. In this treatment, the surrounding fluid 
was replaced by random forces, which only act on the surface of clusters. This situation 
is equivalent to a Brownian motion described by the stochastic Langevin equation, where 
random forces replace the collisions with the fluid, and a constant viscous force represents 
the drag of the cluster motion immersed in the fluid. These two types of forces are related 
by the fluctuation-dissipation theorem. Physically, the interaction between the clusters and 
surrounding liquid is not fully stochastic, the random forces using in this method should be 
carefully. 

For the same purpose, Gong and Sun proposed an alternative CPMD for finite system.- 
In their approach, the system Lagrangian is extended to include the PV term, where P is 
the external pressure and V is the volume of nanoclusters. By writing the volume as a 
function of atomic coordinates, the constant pressure can be readily achieved without any 
pressure-transmitting liquid and without any artificial parameters. Hereafter this method 
is named as extended method. In the application level, the key issue of this method is to 
express the cluster volume as a proper function of atomic positions. Since without periodic 
boundary conditions, even without a regular shape, it is nontrivial to get a proper definition 
of volumes. Gong and his co-workers have proposed a few definitions. In the original CPMD 
paper, 9 Gong et al decomposed the cluster volume as the summation of individual atoms, 
which has been used for metallic systems. Another definition due to Gong and his co-worker 
is the cluster approximated by an ellipse, thus the volume can be written in term of the 
principal radii of gyration. To calculate the enthalpy of clusters, Calvo and Doye^ give a 
more precise definition of volume as the minimum polyhedron enclosing the cluster. 

Carefully using these new methods, now it is possible to make theoretical calculation 
for the low-dimensional system under pressure.- Over the past years, the extended method 
has been used to study the structure and elastic properties of silicon clusters,— metallic 
clusters,— nanotubes,— CdSe nano crystals,— 6*60 and diamond clusters,— etc. The di- 
rected method has been successfully used for the nano systems including the silicon clusters,- 
nanotube,— CdSe nanocrystal,— Cqq and diamond cluster,— etc. The Stochastic method has 
been employed to the Au clusters. 8 

In this paper, we have reviewed the current understanding about the extended method, as 
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well as several recent theoretical advances of nano systems under pressures. The technique 
details of the extend CPMD (ECPMD) are presented in section II; The numerical tests of 
ECPMD are presented in section III; In section IV, the definition of the volume for nano 
systems are recalled; In section V, the application of the new method for some nanocrystals 
and nanotubes are presented; Finally, we summarize the major conclusions in section VI. 

II. The extended constant-pressure MD (ECPMD) method 

Considering a real iV-atom system, its Lagrangian L rea i takes: 

N D 2 

L real = £ ^ " #{*}) (1) 

i 1 

where r\, rrii, Pi are the position, mass, momentum of ith atom respectively, and <fi is the 
interaction potential. In the extended method introduced by Gong and Sun,- the system is 
extended to include a PV term, and Lag rangian L ex tend reads: 

N D 2 

L ex tend = £ ^~ ~ WW) + P^V) (2) 

i Arm 

where V and P ext are the volume of the system and the external pressure respectively. 

The equations of motion(EOM) for the extended system derived from the Lagrangian 
Lextend are, 

d / 9 Lextend \ QL extend /q\ 

The forces acting on the atoms compose of two parts, i.e., the force due to interatomic 
potential {~T j) and the one due to the PV term ( / PV ). EOM derived from Eq. 3 produces 
the constant pressure ensemble for the real systems, which can be obtained according to 
virial theorem. 

I N N N 

< ^7(£ m i v i ~ E r i • - £ r« • PextViV) >= (4) 
where Vi is the velocity of ith atom, and <> denotes average. Then we have, 

Af Af Af 

< £ rmvi - £ ri ■ >=< £ r i ■ PextS/iV > (5) 

i i i 

In the classic statistical physics, there is a basic assumption that any statistical result 
can be obtained exactly from Newton's mechanics, and any statistical quantity should be a 



function of coordinates and velocity of atoms.— Obviously the volume can be written as a 
cubic homogeneous function of atomic positions. 

V = V( ,rf, ). (6) 

where is the position of the ith atom. In fact, all the occupied space by a cluster can 
be divided up into tetrahedra with atoms at their corners, thus the volume of a cluster is 
the cubic homogeneous function of three Cartesian components of atomic positions. Let 
(xi, yi, Zi)(i = ~ 3) denote the Cartesian coordinates of the four vertexes of a tetrahedron, 
thus the volume can be written as, 



^tetrahedron „ 




(7) 



xo -x x y - yi z - z\ 
x -x 2 2/o - 2/2 z - z 2 
xo -x 3 y - 2/3 z - z 3 
The total volume is the summation of the each individual tetrahedron, V = V te trahedron- 
According to Euler theorem, 

$>- v ^ = 3y (8) 

i 

Finally, we end up with 

I N N 

Pext = Pint =< 7777 (E m i V i ~ E r i ■ V*0) > (9) 
i i 

where P{ n t refers to the internal pressure, since the external pressure P ex t is a constant, P{ n t 
is also a constant. Thus, by writing the volume as a function of atomic coordinates, the 
constant-pressure MD is achieved. 

In some special cases, the system size in certain directions fixed( in other words, the 
volume of systems is independent of the atomic position along the special direction). For 
this system, external pressure corresponds to an uniaxial pressure. The uniaxial pressure 
can be realized by only including the one or two components of ~J" PV . For example, If only 
x-component of / PV includes in the simulation, this means the volume is independent on 
the y and z components. Now equation 5 becomes 

N N N Qy 

< ™>iVli - YI X * ' V»0 > = < Y X iPextQ— > (10) 

i i i 1 

Then we have, 

I N N 

Pext = Pxx =< y(Y m i V xi ~ Y T i ' V«0) > (H) 
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where P xx is pressure along x-direction. 

Similarly if only x and y-component of / PV includes in the simulation, this means the 
volume is independent on the z component. Now we have, 

I I N N 

Pext = APxx + Pyy) =< nf^E m i^li + ~ X i ' V«0 + Yi ' V*0) > (12) 

i i 

The above equation has been used in the study of nanotubes under radial pressures.— 
It also can be used for the surface systems. 

Combining with the constant temperature method, the constant pressure method could be 
readily extended to constant-temperature and constant-pressure ensemble. This extension 
is straight forward. Most simulations at finite temperature in this paper are preformed by 
combining with a Nose-Hoover thermostats.— The extension of the method to ab — initio 
molecular dynamics is also simple.— 

The ECPMD method is different from Andersen-Parrinello- Rahman CPMD(APR- 
CPMD) physically, which has been misunderstood by a few authors. First all, the volumes 
in ECPMD and APR-CPMD play the different role. In APR-CPMD method, the volume is 
a generalized coordinates, which has equal importance as an atomic coordinate. However in 
ECPMD scheme, the volume is just a function of atomic coordinations, even not a dynamics 
variable. Secondly, the Lagrangian includes a virtual kinetic energy and mass associated 
with the volume in APR-CPMD, which is absented in ECPMD, thus the atomic dynamics 
in the methods could be different. Finally, in APR-CPMD, the responding of system to 
the external pressure is essentially linear and global, i.e., all the atomic position is linearly 
scaled in the same time. However, in ECPMD, the responding of systems to the external 
pressure is truly local and non-linear. This is especially important for the inhomogeneous 
system. It needs to point that, although the two approaches could be different in dynamics 
level, in the thermodynamics level, both do realize the constant pressure ensemble. 

III. The volume of a cluster. 

One of key issues of the ECPMD is to properly define the volume for a finite system. 
To do this, one should keep two points in mind. One is the intrinsic uncertainty due to 
thickness of cluster surface, which depends on what kinds of materials used to explore its 
thickness. The previous studies on the cluster—, nanotube and nanowall^ have met this 
problem. Another one comes from the computational consideration. Geometrically, one can 
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calculate the volume for any cluster, but it is non-trivial to find one easy to implement and 
computational cheap. 

Before discussing the specific definition of volumes, we would like to make some general 
comments. First, the force due to the PV term just acts on the surface atoms, because 
only the motion of surface atoms directly changes the volume of systems. This is consistent 
with the fact that the pressure-transmitting liquid only interacts with the surface atoms. 
Secondly, the calculation of volume could be much different specific forms, however as only 
as each different form gives the same volume for all the configurations, these forms will 
produce the same dynamics. This is easy to understand mathematically. 

More generally, let V\ and V 2 to be two different definitions, and V 2 —aV\ + b, where a and 
b are constants. The partition function calculated by V\ and V 2 has following relationship, 

Z 2 (P) = J e-n*+ PV » X\dn = J e -/W(«v 1+ &)) jj rfr . 

= J e -n*+*PiV 1+ b/a)) -Q dn = e ~(3Pb/a J e - W+ aPVl) J] fo. = e'^Z^aP) 

Supposing Ai and A 2 are the ensemble average of a physical quantity obtained by using 
V\ and V 2 respectively. According to statistical physics, the thermal average of a physics 
quantity A reads, 

MP) = yJ M{ri})e-^ +PV ^ Rdn, (13) 

A 2 (P) = ±-J A{{n})e~^ PV ^\{dn = ±-J A({r t })e^ +a ^n^ (14) 

where /3=1//cbT, Z\ and Z 2 are the partition function corresponding to V\ and V 2 respec- 
tively. Comparing above two equations, one can easily conclude that 

A 2 (P)=A 1 (aP). (15) 

It implies the physics could be the same for the two different definition of volumes, but it 
may happen in different pressures, if the two definitions have the linear relationship. Eq.15 
also provides a very useful tool for comparison MD results, where different volumes are used. 

If the volume uncertainty due to surface could be neglected, the exact volume can be 
calculated in principle. One of the very accurate definition of volume is writing the volume 
as the summation of all the no cross tetrahedron formed by four atoms, 

V — ^ Vtetrahedron, (16) 
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where Vtetrahedron can be calculated using Eq.7. Although, this scheme for volumes is exact, 
in computational view, it brings large overloading for simulations. In fact, this method was 
not found in any real simulations. 

From computational viewpoints, one usually needs to find a more reliable and cheap way 
to calculate volumes. One of the simple and sufficient ways is to approximate the volume of 
each atom based on the Wigner-Seitz sphere, i.e, the scaled volume of the atomic sphere to 
replace the Wigner-Seitz primitive cell (hereafter labeled as VWS), which has the following 
form, 



here r c keeps between 1st and 2nd nearest neighbors, iVj is the numbers of the nearest 
neighbor of the ith atom, and the summation runs over all the first nearest neighbors of the 
ith atom, ji is a scale factor. For close packed structure, ji is approximated to have the 
value of 1.353. 

VWS was found to work well for metals. Fig.l (middle panel) shows the exact volume and 
one calculated by VWS for bulk Ni liquid at 3000K and 5GPa. One can see that the volume 
calculated by VWS does not recover the exact one instantly. However the instant fluctuation 
can be much reduced by a short time average. We find that the short time average is in 
excellent agreement with the exact one (up panel of Fig.l). Since most physical quantities 
are calculated through time average, we believe the instant fluctuation will result in little 
effect on physical results. 

For most clusters, the ellipsoid is a good approximation to its shapes, its volume can be 
also approximated by the volume of the ellipsoid. The volume of an ellipsoid is determined 
by three semi-axes, which can be given by the radii of gyration Ri(i = 1,2,3) of this cluster. 
This definition of volumes was first used in studying the glass transition of Al n clusters by 
Sun and Gong,— and recently extended by Baltazar et al..— According this definition, the 
volume of the cluster is, 



Where C is a scaling constant, which can be adjusted appropriately according to its real 
volume. Following Baltazar et al.,— the volume can be re-expressed as, 




(17) 



V = C—R1R2R3 



(18) 




(19) 
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where det(J) is the determinant of the inertia tensor /. This definition was found to work 
quite well for Cqo and Si nanocrystal.— It is also recommended for metal systems. 

As we mentioned above, the volume is only determined by the position of surface atoms, 
thus the volume can be written as the minimal polyhedron enclosing the finite system. This 
definition have been used by Calvo and Doye to calculate the cluster enthalpy 1 - In using this 
method, one should pay special attention for clusters with large negative surface curvature, 
since it is easy to judge the surface atoms for a cluster with positive structural curvature, but 
it may be much subtle for part of clusters with negative structural curvature. For finding the 
minimal polyhedron, the most used one is called the quick convex hull algorithm.— Recently, 
this approach has been used for studying the structure transition of CdSe nanocrystal.— 

For some special structures, a specific definition of volume will much simply the comput- 
ing. For example, the surface atoms of nanotubes and fullences can be easily located, the 
volume of a nanotube can be defined through the minimal polyhedron method. 

IV. The numerical tests of ECPMD 

The equation 3 does produce the constant pressure ensemble as shown below. As an 
example, we simulated the carbon nanotube(CNT) and Cqq at 300K based on equation 
3. In this study, the volume is defined through the minimal polyhedron method, and the 
temperature is maintained by Nose-Hoover thermostate.™ The interaction between carbon 
atoms is described by a parameterized many-body potential.—^ In the calculation for CNT, 
the pressure is only applied to all directions normal to axis. In Fig. 2, we present the volume, 
enthalpy and pressure as a function of times for Cqq. From this figure, we can see that the 
evolution of the instantaneous volume, pressure and enthalpy fluctuates around the average 
value, and the average pressure equals to the applied external pressure. The correlation 
between the volume and pressure can also be clearly observed. The similar results are 
shown in Fig. 3 for CNT. The equation 3 now has been tested in many finite systems, and 
in all cases, the constant pressure ensemble is guaranteed. 

To show that the external pressure equals to the internal pressure (Eq.9 and 12), Fig. 4 
shows the internal pressure as a function of external pressure for both Cqq and CNT. For 
C*6o and CNT, the internal pressure is defined as Eq.9 and 12 respectively. The simulation 
results clearly show that the constant pressure for both cases hold. For other system, the 
similar results have been obtained. 
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V. The applications of ECPMD 



Materials under pressure have plenties of phenomena and attract people for hundreds 
of years. The high pressure experiments provide very important information relevant to 
the structure stability and bonding of materials. Recently the studies on nano-systems un- 
der pressure show fruitful new phenomena.— i22j22j2L22j22jM Promoted by the high pressure 
experimental work, Gong and his coworkers have studied the structure and properties of 
nanosystems under pressures by using ECPMD. In the following of this review, the theo- 
retical approaches for the finite systems under pressure based on the ECPMD are recalled, 
which includes, 

(a) The elastic properties and melting behavior of metallic nanocrystals 

(b) Structure transformation of CdSe nanocrystals 

(c) Pressure induced hard-soft transition of carbon nanotubes 

(a) The elastic properties and melting behavior of metallic nanocrystals 

In the studies for metallic systems, the well-tested many-body potentials are used, namely 
glue potentials for Au,— Sutton-Chen potential for Ni,— and the tight-binding model for 
Ag.— The volume is calculated basing on VWjS', the system temperature is realized by using 
Nose-Hoover thermostats. 

The bulk modulus is one of the most important parameters of materials, which reflects 
the the elastic properties of materials. The elastic properties of Au, Ag and Ni nanocrystals 
have been studied by using ECPMD. Fig. 5 shows the pressure and energy as a function 
of reduced volumes for Ni nanocrystals at 300K, and the counterpart of the bulk phase 
calculated by APR-CPMD. In this figure, the energy is relative to the minimum energies, 
and the volume is renormalized by the equilibrium volume at zero pressure and 300 K. 
Clearly nanocrystals is softer than the bulk phase, which is reflected by the larger volume 
change for nanocrystal than bulk for the same applied pressure. The similar results are also 
found for other metallic nanocrystals. The bulk modulus can be obtained by fitting the 
energy- volume plot. The obtained bulk modulus as a function of the size of nanocrystals is 
shown in Fig. 6, where the solid line is the linear fitting. It can be seen that the MD data 
follows a straight line quite well, which implys the elastic constants are reduced inversely 
with the size of nanocrystal, similar to many other properties for nanosystems. The bulk 
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modulus for different temperature is also obtained by fitting to the energy- volume plot. The 
obtained bulk modulus as a function of temperature is shown in Fig. 7. The present results 
show that, the bulk modulus decreases with the increase of temperatures, which is similar 
to the bulk phase, and also consists with the basic thermodynamics results. 

The melting behavior is one of the common phenomena in nature, which is also one of the 
most important process relevant to the properties of materials. The basic thermodynamics 
shows that the melting temperatures are strong affected by the external pressures, which is 
characterize by the so-called Clapeyron equation for bulk materials. Although the melting 
behavior of nanoclusters have been wildly studied over the past decades, the pressure effect 
on the melting behavior was not well understood for nano systems yet. Recently Ye et al. 
have carried out a detailed study for the melting behavior of Ni nanoclusters under pressures. 

Fig. 8 shows the melting points(T^) as a function of pressures for Ni 5G i. Consistent 
with the basic thermodynamics, the melting temperatures is increasing with the increasing 
of pressures. The similar results is found for Nii^ and bulk materials. The latent heat 
versus pressure for Ni^i cluster is shown in Fig. 9. From this figure, we can see that the 
latent heat seems to be a constant in the studied range of pressures, where the average 
value is about 0.0806eV. The volume difference between the solid state and liquid state at 
the melting point versus pressure for Ni^i cluster shown in Fig. 10. Assuming the latent 
heat is the constant, The melting temperature and volume difference are related through 
the Clapeyron equation quite well. 

(b) Structure transformation of CdSe nanocrystals 

Ye et al studied the structure transformation of CdSe nanocrystals using the ECPMD. 
The empirical potential developed by Rabani.^A has been used to describe the interatomic 
interaction. Most of their simulations are carried out at 300 K by using a Nose-Hoover 
thermostat.— The remarkable structure character of nanoclusters is the large surface- volume 
ratio, thus it can be expected that the surface could play an important role for the structure 
transformation in nano systems. In order to study the effect of the surface structure on 
the transition mechanism, they use nanocrystals of two different shapes, i.e, the spherical 
and faceted one, consisting of 500 to 5000 atoms. The initial configuration of the spherical 
nanocrystal is simply cut from the bulk CdSe of WZ structure. Faceted nanocrystals with 
well-defined surface structure are obtained by cleaving the bulk lattice along equivalent 
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(100) WZ planes and at (001) and (OOl)planes perpendicular to the [001] direction of the c 
axis. The volume of the nanocrystals is approximated based on finding the subset of atoms 
forming the smallest convex polyhedron. 

Ye et al have observed the transformation from wurtzite to rocksalt structure, but the 
process of transformation is strongly dependent on the shape and size of the nanocrystals. 
Upon loading the pressure, the spherical CdSe nanocrystals is found to directly transfer 
to rocksalt structures with nanoscale grain boundary formed, while the faceted ones can 
first transfer to hexagonal MgO structure, and then the final rocksalt structure with grain 
boundary free. These results are similar to that calculated by the direct method for the 
same systems. 

Fig. 11 shows their calculated the volume-pressure plot. From this figure, it clearly in- 
dicates the structural transformation of the CdSe nanocrystal up loading pressures. The 
volume of the spherical Cd 502 Se 50 2 nanocrystal decreases smoothly with increasing pressure 
up to a critical pressure~ 8.0 GPa, at which the volume decreases abruptly as a result of 
the transformation from WZ to RS(left penal of Fig. 14). This is in good agreement with 
the high pressure experiment for the same system.— For faceted nanocrystals, an intermedi- 
ate structure(C— >D) clearly exists between WZ and RS. Detailed analysis of the variations 
of coordinations shows that the WZ structure of faceted nanocrystal transforms to a five- 
fold coordinated structure around 1.4GPa. The five-fold coordinated structure has been 
reported as a stable phase of MgO under hydrostatic tensile loading. 37 ' 38 When the pressure 
continues to increase, the five-fold coordinated structure transforms to six-fold coordinated 
RS structure. The transition process is also found to be highly hysteretic. Upon pressure 
releasing, The rock salt structure remains stable down to pressures significantly below the 
observed "upstroke" transition pressure. As low as 0.5 GPa, the sample begin to restore 
WZ structure. At atmospheric pressure, the WZ is recovered but with a few defects near 
surfaces. 

For spherical nanocrystals, the transformation pressure decreases with size increasing. 
This trend coincides with the experimental results of Alivisatos.— ^ (Note, in Ye et a/'s 
studies, "upstroke" transition pressure is used as the transition pressure) In contrast to the 
spherical one, the transformation pressures (both from four-fold to five- fold and from five- 
fold to six-fold one) of facet nanocrystals increase with the increasing of crystal sizes. Fig. 
12 and Fig. 13 shows the transformation pressure as a function of sizes for facet nanocrystals. 
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The different dependence of transformation pressure on the size has been discussed based 
on thermodynamics considerations by including the surface effect. 

For all spherical nanocrystals, nano-scale 'grain' boundaries are formed during and after 
the transformation(see Fig. 14). As the size of the nanocrystal increases, the multiple grains 
phenomena become more and more obvious. In contrast, the facet one is almost grain 
boundary free. The generation of grain boundary has been discussed based on the nucleation 
mechanism.— 

(c) Pressure induced hard-soft transition of carbon nanotubes 

Gong and his co-workers have presented a detailed investigation on the behavior of car- 
bon nanotubes under hydrostatic pressures by ECPMD.— d&M In their simulation, a few 
single- walled carbon nanotube(SWCNT), double- walled carbon nanotubes(DWCNT) and 
multi-walled carbon nanotube(MWCNT) are investigated. For DWCNT, both commen- 
surate and incommensurate one are considered. DWCNTs consisting of tubes, which has 
the same chirality, are typically commensurate; otherwise they are incommensurate. The 
volume is calculated by the minimal polyhedron method. The periodic boundary condition 
in the axial direction, and free boundary condition in the radial directions are used. The 
interaction between carbon atoms is described by a parameterized many-body potential.— 1 ^ 
The intertube and intratube van der Waals interaction are modelled by the Lennard- Jones 
(LJ) potential.— ^2 To confirm the results of the classical molecular dynamics method, they 
also have repeated some calculations by ab-initio molecular dynamics method. 

They found that all studied nanotubes(SWCNTs, DWCNTs and MWCNTs) undergo a 
pressure- induced hard-to-soft phase transition. The hard phase at low pressure exhibits a 
typical bulk modulus of 100 GPa, while the soft phase at high pressure exhibits a bulk 
modulus of only ~1 GPa. Fig. 15 shows the pressure and the total energy as a function of 
reduced volume for a (10,10) nanotube at 300 K, where the energy at zero pressure is set to 
zero. Clearly, a transition at ~1.0 GPa is observed. Below the transition, the hard phase 
has a radial compressibility of 0.01 GPa -1 . Above the transition, the soft phase has a radial 
compressibility about two orders of magnitude larger. The similar behavior was observed 
for other tubes. 

After the hard-to-soft transition, the cross section of nanotubes changes from circular 
to elliptical shape. The evolutions of the cross-section shape, the bond length, and the 
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bond angle with increasing pressure for a (10,10) nanotube are shown in Fig. 16. where 
two principal axes (long axis a and short axis b) are used to characterize cross section. 
Below the transition pressure, a remains almost equal to b, defining a circular shape( see 
Fig. 16). Above the transition pressure, a becomes larger than b, defining an elliptical shape. 
Eventually, as the long axis continues to increase and the short axis continues to decrease, 
the elliptical shape undergoes another transition to a dumbbell shape. Under even higher 
pressure, the dumbbell tube can become so flat that the spacing between the opposite side 
walls approaches the layer spacing in the graphite (~3.35A). 

The trend of change in bond length and bond angle provides a good explanation of the 
hard-to-soft transition. Fig. 16(b) shows that the percentage change in bond length and 
bond angle increases simultaneously with increasing pressure below the transition, indicating 
a uniform shrinking of the circular shape under pressure. Above the transition, the bond 
length remains unchanged but the change of bond angle increases sharply with increasing 
pressure. Since it costs much more energy to change bond length than to change bond 
angle. Below the transition, the structural response to the external pressure is largely taken 
by the changing bond length of a circular shape, giving rise to a hard phase; while above the 
transition, the structural response to the external pressure is largely taken by the changing 
bond angle of an elliptical shape, giving rise to a soft phase. 

The critical transition pressure depends strongly on the tube radius. Fig. 17 shows the 
simulated transition pressures as a function of tube radius for SWNTs(solid dots). The 
smaller the radius, the higher the transition pressure. To understand the above simulation 
results, they also provide a general analysis based on continuum elastic theory. According 
to their deduction, transition pressure P t has, 



where D is the constant related to the elastic properties of NTs, and R is the tube radius. 
This analytical dependence of Pt on Rq is in very good agreement with the MD simulations, 
as shown in Fig. 17. 

The bulk modulus of the hard phase follows, 



Pt 



3D 



(20) 




C 



(21) 



2R 
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where C is another constant related to the elastic properties of NTs. This analytical de- 
pendence of bulk modulus on Rq are in very good agreement with the MD simulations, as 
shown in Fig. 17. 

The similar pressure-induced structural transition has been found for all studied DWNT's. 
Comparing with SWNT, the transition pressure of DWNT is much enhanced, but it still 
follows the Eq. 20. Figure 18 presents the transition pressure as a function of the both inner 
and outer tube radius. The results imply that the van der Waals forces between two tubes 
does affect the transition pressure. The transition pressure of a outer tube in DWNT can 
be increased largely by inserting an inner tube. In fact, the DWNT can be considered as a 
psudo-single-walled nanotube with effective thickness, which should be larger than the real 
SWNT. 

The remarkable feature of MWNTs is its encapsulation effects, especially, when the sys- 
tem is undergoing pressure, the outer shell acts as protector for the inner shell. They found 
that, the response pressure of inner tube is much smaller than the external pressure, while 
the response pressure of outer tube is much closer to the external pressure. To characterize 
the pressure transmission, Ye et al define a response pressure for the tubes.— Left panel 
of Fig. 18 shows the response pressure as functions of external pressure of (5,5)@(10,10) 
DWCNT. Form this figure, on can see that the response pressure of both inner and outer 
tube increases linearly with the external pressure below 8 GPa, while the value for inner tube 
is about two times smaller than outer one. Interestingly, for all the studied tubes, when the 
external pressure is higher than a certain value, at which the hard-to-soft transition happens, 
(~8 GPa for (5,5)@(10,10) and ~7 GPa for (5,5)@(10,10)@(15,15)), the response pressure 
of inner tube increases sharply, while increasing of the response pressure of the outer tube 
slows down. 

Ye et al assume a linear relationship between the response pressure and the external 
pressure before the structural transition happens, they define pressure transmission efficiency 

(3 by 

p r = a + f3P e 

where P r and P e are the response pressure of the inner tube and the external pressure 
respectively, a is the response pressure of the inner tube without external pressure. Fig. 19 
presents the pressure transmission efficiency of (n,n)@(n+5,n+5) DWCNTs as a function of 
the radius of the outer tube. The pressure transmission efficiency is found to increase with 
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the tube radius. 

In contrast to commensurate DWCNTs, the incommensurate DWCNTs have lower trans- 
mission efficiency. The pressure transmission efficiency are 0.30 and 0.35 for (6,6)@(19,0) and 
(10,10)@(26,0) respectively, while the pressure transmission efficiency of their commensurate 
counterpart (6,6)@(11,11) and (10,10)@(15,15) are 0.35 and 0.43 respectively. Obviously the 
morphology combination does affect the vdW interaction between inner and outer tubes. 
The calculations show that the pressure transmission of the commensurate DWCNTs is 
more efficient. This might be due to the fact that the atomic positions in adjacent shells 
are well matched in commensurate DWCNTs, meanwhile the intralayer vdW force favors 
commensurate tubes.— 

VI. Summary 

By writing the volume of a system as a function of coordination of atoms and extending 
the Lagrangian of the system to include a PV term, a constant-pressure molecular-dynamics 
method can be achieved in a simple but physically rigid way. This method is different the 
traditional constant-pressure one by treating volume as a part of potential in steady of a 
generalized dynamics variable. This method is specially suitable for finite systems and the 
system without periodic boundary conditions. In this paper, the varies of application and 
some technique key issues of this method are reviewed. The method is fairly general and 
can find widespread applications. 
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FIG. 1: Bottom panel: The volume per atoms for liquid Ni calculated by VWS (thin line) and the 
exact one (thick line) obtained by traditional constant pressure molecular dynamics simulation at 
3000K and 5GPa. Up panel: the same as bottom panel except the volume averaged over a short 
time. The two lines are close to each other 
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FIG. 2: The evolution of the instantaneous internal pressure (middle penal), enthalpy (low penal) 
and volume (up penal) during ECPMD runs for Cqq. The ECPMD does recover a constant pressure 
simulation. 
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FIG. 3: The same as Fig. 2 except the system is carbon nanotube. 
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FIG. 4: The internal pressure as a function of external pressure for C$o and carbon nanotube. 
Clearly they are equal! 
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FIG. 5: The equations of state for Ni nanocrystals (filled circles) and bulk phase (open circles), 
where the data for iV^isi nanocrystals and bulk phase are calculated by ECPMD traditional 
molecular dynamics simulation respectively. 
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FIG. 6: The bulk modulus as a function of the size of nanocrystals. Circles: calculated data, line: 
the linear fitting to the data. 
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FIG. 7: The bulk modulus as a function of temperature for iViiggs. 
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FIG. 8: The melting points(T^) as a function of pressures for Ni^Qj. The melting temperatures 
is found to increase with the increasing of pressures, in agreement with the basic thermodynamics. 




FIG. 9: The latent heat versus pressure for Ni^i cluster. The latent heat seems to be a constant 
in the studied range of pressures. 
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FIG. 10: The volume difference between the solid state and liquid state at the melting point 
versus pressure for Ni$qi cluster. Dash line: deduced from Clapeyron equation, Square: from the 
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FIG. 11: Volume versus pressure for faceted Cdu^Seu^ (right) and spherical Cd^Se^ 
nanocrystal (left). The discontinuity of the slopes indicates structural transformation. 
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FIG. 12: Variation of the transformation pressure with radius for spherical CdSe nanocrystals at 
300 K. With increasing nanocrystal size, the transformation pressure decreases. 
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FIG. 13: The transformation pressure as a function of radius for facted nanocrystals at 300 K. In 
contract with spherical one, the transformation pressure increases with increasing nanocrystal size. 
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FIG. 14: Domains after the structural transformation in spherical nanocrystals. (a), (b), (c) and 
(d) show the cross sections through the middle of the nanocrystals of radius 19 A, 23 A, 27 A and 
33 A, respectively. The grain boundaries are shown with gray atoms. 



29 



o 

S- 

s 

S- 




-2.0 




0.80 0.90 1.00 
Reduced Volume 



1.10 



FIG. 15: The energy and pressure as a function of the reduced volume for (10,10) carbon nanotube 
at 300 K. The minimum energy is set to zero, and the volume is normalized by the equilibrium 
volume without the external pressure. At about 1.0 GPa, the hard phase with bulk modulus of 
about 100 GPa transforms into the soft phase with bulk modulus of just a few GPa. 
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FIG. 16: The length of the long and short function of pressure for (10,10) nanotube. The 

shape of cross section at some selected pressures is plotted at the bottom of the figure. The absolute 
relative change of bond length and bond angle as a function of pressure for (10,10) nanotube, the 
data is obtained by quenching the system from 300 K to K at constant pressure. 
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FIG. 17: The transition pressure (upper panel) and the elastic modulus (lower panel) as a function 
of tube radius at 300 K. The solid line is a least-square fit to the data using Eq. 20 (upper panel) 
and Eq. 21 (lower panel). The simulated data follows nicely with the predicted behavior. 
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FIG. 18: Transition pressure as a function of radius of nanotubes at 300 K for a few DWCNTs. 
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FIG. 19: Response pressure as a function of the external pressure for (5,5)@(10,10) DWCNT(left) 
and (5,5)@(10,10)@(15,15) TWCNT(right). The system response pressure is exactly the same as 
the external pressure. The response pressure of inner tube is much smaller than that of the outer 
one before structural transition occurs. 
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FIG. 20: Pressure transmission efficiency of commensurate (n,n)@(n+5,n+5) DWCNT (with 
n=5,6,7,8,9,10) versus the outer tube radius. The transmission efficiency increases with the tube 
radius. The incommensurate DWCNTs (6,6)@(19,0) and (10,10)@(26,0) is also shown for compar- 
ison, of which the transmission efficiency is smaller. 
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